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Abstract 



We prove a Korn-type inequality in H(Curl; SI, R 3x3 ) for tensor fields P mapping n to R 3x3 . More 
precisely, let SI C R 3 be a bounded domain with connected Lipschitz boundary dSl. Then, there exists a 
C*j ' constant c > such that 

3 : 

1—5 . c I p lL2cn,a3x3) < ||symP|| L2(fJiR 3x 3) + |CurlP| L2{niE 3x 3) (0.1) 

o . 

. holds for all tensor fields P 6 H(Curl; SI, R 3x3 ), i.e., all P £ H(Curl; Si, R 3x3 ) with vanishing tangential trace 

on dSl. Here, rotation and tangential trace are defined row-wise. For compatible P, i.e., P = Vti and thus 
CurlP = 0, where v G H^^R 3 ) are vector fields having components v n , for which Vv n are normal at dSl, 
the presented estimate (|0.1|) reduces to a non-standard variant of Korn's first inequality, i.e., 



Oh 

< 



o 



X 



c l|Vw|| U 2 (niK 3x3) < |symVu| L 2(n iK 3xs) 



On the other hand, for skew-symmetric P, i.e., symP = 0, 1)0. l|l reduces to a non-standard version of 
Poincare's estimate. Therefore, since (|0. 1 f) admits the classical boundary conditions our result is a common 
generalization of the two classical estimates, namely Poincare's resp. Korn's first inequality. 
Key Words Korn's inequality, gradient plasticity, theory of Maxwell's equations, Helmholtz decomposition, 
s ! ' Poincare/Friedrichs type estimate 

; 1 Introduction: Infinitesimal Gradient Plasticity 



The motivation for our new estimate is a formulation of infinitesimal gradient plasticity [5] • Our model is taken 
from Neff et al. [9]. Let !1 C I 3 be a bounded domain. The goal is to find the displacement u : [0, oo) x SI h-> R 3 
and the possibly non-symmetric plastic distortion tensor P : [0,oo) x f2 H> R 3x3 , such that in [0,oo) x SI 



Div a = f, a = 2pi sym(Vu - P) + A tr(Vu - P) id, 

P £$(£), S = a - 2/j, sym P - Curl CurlP, (1.1) 

hold. The system is completed by the boundary conditions 

w(*,a;)=0, v{x) x P(t,x) = V (t, x) G [0, oo) x OS} 

and the initial condition P(0,x) = for all x G SI. The underlying thermodynamic potential including the 
plastic gradients in form of the dislocation density tensor Curl P is 



I 

-In 



fi\ sym(Vu - P)| 2 + ^| tr(Vu - P)| 2 -f-u + fi\ symP| 2 + ^L 2 C \ CurlPp 



Here, fi, A are the elastic Lame moduli and a is the symmetric Cauchy stress tensor. The system is driven 
by nonzero body forces denoted by /. The exterior normal to the boundary dSl is denoted by v and the plastic 
distortion P is required to satisfy row-wise the homogeneous tangential boundary condition which means that 
the boundary dSl is a perfect conductor regarding the plastic distortion. Q 

Moreover, $ : R 3x3 H> R 3x3 is the monotone, multivalued flow-function with $(0) =0 and $(R 3x n 3 ) c R^- 
In general, E is not symmetric even if P is symmetric. Thus, the plastic inhomogeneity is responsible for 



*This homogeneous tangential boundary condition on P is consistent with v X Vu = on 9f! which follows from u = on c*f!. 
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the plastic spin (the possible non-symmetry of P). The mathematically suitable space for symmetric plastic 
distortion P is the classical space H(curl;f2) for each row of P jTTJ [2J. This case appears when choosing 

In the large scale limit L c — > we recover a classical clasto-plasticity model with local kinematic hardening 
and symmetric plastic strain e p := symP, since then P € Rgym- 

Uniqueness of classical solutions for rate- independent and rate-dependent formulations of this model is shown 
in [9]. The more difficult existence question for the rate- independent model in terms of a weak reformulation 
is addressed in [5] . First numerical results for a simplified rate- independent irrotational formulation (no plastic 
spin, i.e., symmetric plastic distortion P) are presented in cf [17j . In [3J the model has been extended 
to rate-independent isotropic hardening based on the concept of a dissipation function defined in terms of 
the equivalent plastic strain. From a modeling point of view, it is strongly preferable to again have only the 
symmetric (rate) part of the plastic distortion P appear in the dissipation potential. 

The existence and uniqueness can be settled by recasting the model as a variational inequality, if it is possible 
to define a bilinear form which is coercive with respect to appropriate spaces. This program has been achieved 
for other variants of the model in (3J. It had to remain basically open for the above system (jl.lj) . In this case, 

o 

the appropriate space for the plastic distortion P is the completion H sym (Curl; fi) of the linear space 

{P G C°°(n,R 3x3 ) : P n normal at dfl, n = 1,2,3} 
with respect to the norm ||| • ||| , where P n are the columns of P T and 

2 (o) + |CurlP|| L2(f 



|||P||| 2 := | Sym P| 2 +|CurlP| 2 2 



Despite first appearance, this quadratic form indeed defines a norm as shown in [5]. Thus H sym (Curl; Q) is a 
Hilbcrt-spacc. However, in this space it is not immediately obvious how to define a linear and bounded tangential 
trace operator. Since only ||symP|| L 2/QS appears, it is also not clear, how to control the skew-symmetric part of 
P. Therefore, the crucial embedding 

H sym (Curl;fi) C L 2 (0) 

is not clear as well. As a consequence of our main result of this paper we obtain that nevertheless 

H syin (Curl;Q) = H(Curl;fi) 

holds with equivalent norms in case the domain Q is simply connected and has a Lipschitz boundary. The result 
of this paper has been announced in [10] , 

For the proof of our main result (|0.ip we combine techniques from electro-magnetic and elastic theory, 
namely the Helmholtz decomposition, the Maxwell compactness property and Korn's inequality Their basic 
variants are well known results which can be found in many books, e.g., in [5] and the literature cited there. 
More sophisticated and related versions are presented, e.g., in [T^J HH E3 HH1 [H] for Maxwell's equations and 
[TJ [5] for Korn's inequality. 

This paper is organized as follows. After this motivation we introduce our notation, definitions and provide 
some background results. In scction|3jwe give the proof for our main estimates. In the last section|4]we establish 
a connection to a related result by Garroni et al. g] for the two-dimensional case. 

2 Definitions and Preliminaries 

Let ft be a bounded domain in K 3 with connected Lipschitz continuous boundary T := dfl. 
2.1 Functions and Vector Fields 

The usual Lebesgue spaces of square integrable functions, vector or tensor fields on Q with values in E, M 3 or 
M 3x3 , respectively, will be denoted by L 2 (SI) . Moreover, we introduce the standard Sobolev spaces 



H(grad;ft) = {u G 1/(0) 
H(curl;0) = {v G L 2 (f2) 
H(div ; r2) = {v G L 2 (ft) 



• gradw G l 2 (fl)}, HlH(Grad;fi) != Hl 2 (Q) + llg rad u ll 2 (Q) 



l(Grad ; n) 

CUrlw G L 2 (fi)}, MlH(curl;Q) != MIl 2 (Q) + ll CUrl V \\l*(f, 

(div;f2) 



divi) G L 2 (fi)}, HlH(div;fi) : = l v l 2 LHn) + ll di H| 2 2(n) 
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H(grad;fi) is often denoted by H 1 (f2). Furthermore, we define their closed subspaces H(grad;f2), H(curl;f2) as 

o 

complctition under the respective norms of the scalar resp. vector valued space C°°(f2) of compactly supported 
and smooth test functions resp. vector fields. In the latter Sobolev spaces the usual homogeneous scalar resp. 
tangential boundary conditions 

u\r — 0, v x v\v = 

are generalized, where v denotes the outer unit normal at T. We note in passing that v x w|r = is equivalent 
to t • w|r = for all tangential directions r at T, which means that v is normal at T. Furthermore, we need the 
spaces of irrotational or solenoidal vector fields 

H(curl ;fi) := {v g H(curl;Q) : cur\v = 0}, 

H(curl ;fi) := {v g H(curl;0) : curlv = 0}, 
H(div ;ft) := {v g H(div;ft) : divu = 0}, 

where the index indicates vanishing curl or div, respectively. All these spaces are Hilbcrt spaces. E.g., in 

o 

classical terms we have v g H (curio; ft), if and only if 

curlu = 0, v X v\r = 0. 

For an introduction of these spaces see [5J p. 11-12, 148] or [5J p. 26]. The most important tool for our analysis 
is the compact embedding 

H(curl; fi) n H(div; Q) ^ L 2 (0), 

which is often referred as 'Maxwell compactness property', see [SI p. 158] and [U [T5J [THl HH HFJ- A first 
immediate consequence is that the space of so called 'harmonic Dirichlct fields' 

H(Q) := H(curl ; Q) n H(div ; fl) 

is finite dimensional. A vector field v belonging to means in classical terms that 

curlw = 0, divw = 0, v x i>|r = 0. 

The dimension of H(fl) equals the second Betti number of fi, see [6j p. 159] and |T3l Theorem 1]. Since we 
assume the boundary T to be connected, there are no Dirichlct fields besides zero, i.e., 

n(n) = {o}. 

This condition on the domain Q resp. its boundary T is satisfied e.g. for a ball or a torus. 

By a usual indirect argument we achieve another immediate consequence, see [SJ p. 158, Theorem 8.9] or J5J 
Lemma 3.4]: 

Lemma 1 (Maxwell Estimate for Vector Fields) There exists a positive constant c m , such that for all vector fields 
v g H(curl;Q) n H(div;ft) 

IMIl 2 (q) < Cm( ll cur MlL2(o) + l divw lL 2 (a) ) 1/2 - 

By definition of the weak divergence, the projection theorem and Rellich's selection theorem [SJ p. 14] wc 
have from p. 148, Theorem 8.3] or [101 Lemma 3.5], [7J Theorem 3.45] 

Lemma 2 (Helmholtz Decomposition for Vector Fields) We have the orthogonal decomposition 

\J(n) = gradH(grad;Q) H(div ;^). 

2.2 Tensor Fields 

We extend our calculus to (3 x 3)-tcnsor (matrix) fields. For vector fields v with components in H(grad; fl) and 
tensor fields P with rows in H(curl; £1) resp. H(div; 0), i.e., 



v n g H(grad; O), P T = [P x P 2 P 3 ], P„ € H(curl; O) resp. H(div; fi) 



4 



Patrizio Ncff, Dirk Pauly, Karl-Josef Witsch 



we define 



Gradv 



grad T ^i 
grad T ii 2 
grad T i> 3 



J v =Vv, CurlP 



curl T Pi' 
curl T P 2 
curl T P, 



DivP 



div Pi 
divP 2 
divP 3 



where J v denotes the Jacobian of v and T the transpose. We note that v and Div P are vector fields, whereas 
P, CurlP and Grady are tensor fields. The corresponding Sobolev spaces will be denoted by H (Grad; SI), 

H(Grad;Sl), H(Curl;fi), H(Curl;fi), H(Curl ;Sl), H(Curl ;Sl), H (Div; SI), H(Div ;Sl). As usual, we denote by 
symP := 1/2(P + P T ) the symmetric part of a tensor P. 

Let us now present our three crucial tools to prove the new estimate. First we have obvious consequences 
from Lemmas [1] and 

o 

Corollary 3 (Maxwell Estimate for Tensor Fields) For all P G H(Curl; SI) n H(Div; SI) 

||P|| L2( o) < Cm{ |CurlP| 2 L2(n) + ||DivP|| 2 L2(n) ) 1/2 . 
Corollary 4 (Helmholtz Decomposition for Tensor Fields) We have the orthogonal decomposition 

L 2 (S1) = Grad H (Grad; SI) © H(Div ; SI). 
The third important tool is Korn's first inequality [6] p. 207] or [TBI p. 54]: 

o 

Lemma 5 (Korn's First Inequality) For all v G H(Grad;Sl) 

||Gradw| L 2 (0) < v/2lsymGradu|| L 2 (n) . 

3 Main Results 

For tensor fields P G H(Curl; SI) we define the semi-norm 

|||P||| :=(||symP|| 2 2w + ||CurlP|| 2 2(n) ) 1/2 . 

O 

Lemma 6 Let c := max{2, \/bc m }. Then, for all P £ H(Curl; Q) 

\p\vm <£\l p \l 

o 

Proof Let P G H(Curl; SI). According to Corollary @] we orthogonally decompose 

P = Grady + Q G Grad H (Grad; SI) © H(Div ; SI). 

o 

Then, CurlP = CurlQ and we observe Q G H(Curl; SI) n H(Div ; SI) since 

Grad H(Grad; fi) c H(Curl ; 0). (3.1) 

By Corollary [3] we have 

\Q\ L * m <c m ||CurlP|| L2(n) . (3.2) 
Then, by Lemma [5] and (|3.2[) we obtain easily 



\\ p \\h(n) = \\ Gl ~& dv + QWl^n) = || Grad w||L 2(n) 



lL 2 (f2) 



< 2 ||symGradt;|| 2 2(0) + ||Q|| 2 2(0) = 2 ||sym(P - Q)|| 2 2(0) + ||Q|| 2 2(0) 

< 4 ||symP|| 2 2(n) + 5 ||Q|| 2 2(n) < 4 ||sym P|| 2 2(0) + 5c 2 n ||CurlP|| 2 2(n) . 

The immediate consequence is 



□ 
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Theorem 7 On H(Curl;f2) the norms \\ ■ In(Curl'f!) 
H (Curl; fi) and 



id II • I are equivalent. In particular, ||| • ||| is a norm on 



3c>0 VPGH(Curl;0) c |P| H(Curl;0) < ||symP|| L2(n) + ||CurlP|| L2 
Setting P := Gr&dv we obtain by Lemma [H] and p. II) 

o 

Remark 8 (Korn's First Inequality: Tangential-Variant) For all v G H(Grad;f2) 

||Gradw|| L 2 (a) < c||symGradt;|| L2(f2) . 



(3.3) 



This is Korn's first inequality from Lemma [5| with a larger constant c. Since the boundary T is connected, 
i.e., H{9) = {0}, we have Grad H (Grad; fi) = H(Curl ;Q). Thus, (O holds for all v G H(Grad;Q) with 

o 

Gradv G H(Curlo;f2), i.e., with Gradi>„ ; n = 1,2,3, normal at T, which then extends Lemma\^ through the 
(apparently) weaker boundary condition. 



4 Two-Dimensions: a Result of Garroni et al. 

Let be a bounded domain in R 2 with connected Lipschitz continuous boundary L, which is equivalent (in K 2 ) 
to the topological property that f2 is simply connected. For tensor fields P : 51 M> R 2x2 we define analogously 
the Curl-operator by 



Curl P = Curl 



Pll P\2 
Pl\ P22 



curl [P n P 12 ] J 
curl [P 21 P 22 f 



di P12 - d 2 P n 
d\ P22 — d 2 P21 



where now curl denotes the two dimensional scalar rotation and CurlP is a vector. With the appropriate 
changes, Lemma [5] and Theorem [7] hold as well. In particular, there exists a positive constant c, such that 



c||P|| L 2 (n) < ||symP|| LS 



(n) 



|CurlP| L! 



(n) 



holds for all P G H (Curl; fl). 

During the preparation of our paper we got aware that a two-dimensional related result may be inferred 
from Garroni et al. @J- Instead of tangential boundary conditions v x P|r = they impose the normalization 
condition 

(4.1) 



skew P = 0. 

Let us define the total variation measure of the distribution CurlP for P e L x (0) by 
|CurlP|o:= sup (P, CoGradw), 2 <^ , CoGradu 



sup (P, CoGradv) L 2( ) 

t)6C 1 (J2) 
i"I L ~ (c2) <l 



d 2 vi - di vi 
d 2 v 2 - di v 2 



We note 



(P,CoGradw) Li 



(n) 



P11 d 2 vi - P12 9i vx + P 2 i 9 2 «2 - P22 di v 2 . 



Using partial integration, i.e., (P, CoGradv) L 2 ( m = (CurlP, 



'L 2 (fi) 



for v E C x (57), it is easy to see that 

|CurlP|n = ||Curl P| L i(Q) if CurlP G L x (0). In [H Theorem 9] it is shown that for 51 having a Lipschitz 
boundary and a special 'slicing' property, there exists a constant c > 0, such that 



\P\vm) < |symP| L2 



(O) 



Curl P| 



holds for all P G L 1 (O) with (|4.1[) . Their proof uses essentially that in R 2 the operators curl and div can be 
exchanged by the simple transformation, i.e., curl[vi,i>2] = div [— v 2 ,vi] . Thus, such a strong result may 
not be true in higher space dimensions N > 3 and it is open whether the normalization condition (|4.ip can be 
exchanged with the more natural tangential boundary conditions. 
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